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Abstract: 

To match the precision of present and future measurements of iy-boson produc- 
tion at hadron colliders, electroweak radiative corrections must be included in the 
theory predictions. In this paper we consider their effect on the transverse momen- 
tum (p T ) distribution of W bosons, with emphasis on large pp- We evaluate the 
full electroweak 0{a) corrections to the process pp — > Wj including virtual and real 
photonic contributions. We also provide compact approximate expressions which 
are valid in the high-energy region, where the electroweak corrections are strongly 
enhanced by logarithms of sjMyj. These expressions include quadratic and single 
logarithms at one loop as well as quartic and triple logarithms at two loops. Numer- 
ical results are presented for proton-proton collisions at 14 TeV. The corrections are 
negative and their size increases with px- At the LHC, where transverse momenta 
of 2 TeV or more can be reached, the one- and two-loop corrections amount up to 
—40% and +10%, respectively. 



March 2007 



1 Introduction 



The study of gauge-boson production has been among the primary goals of hadron 
colliders, starting with the discovery of the W and Z bosons more than two decades 
ago [lj. The investigation of the production dynamics, strictly predicted by the 
electroweak theory, constitutes one of the important tests of the Standard Model. 
Differential distributions of gauge bosons, in rapidity as well as in transverse mo- 
mentum (pt), have always been the subject of theoretical and experimental studies. 

At large pp the final state of the leading-order (LO) process consists of an elec- 
troweak gauge boson plus one recoiling jet. The high center-of-mass energy at the 
Large Hadron Collider (LHC) in combination with the enormous luminosity will 
allow to explore parton-parton scattering up to energies of several TeV and corre- 
spondingly production of gauge bosons with transverse momenta up to 2 TeV or 
even beyond. In this energy range the electroweak corrections from virtual weak- 
boson exchange are strongly enhanced, with the dominant terms in L-loop approxi- 
mation being leading logarithms (LL) of the form a L log 2L (s/Mjy), next-to-leading 
logarithms (NLL) of the form a L log 2L_1 (s/M^), and so on. These corrections, 
also known as electroweak Sudakov logarithms, may well amount to several tens 
of percent [21 El 13 El El E]- A recent survey of the literature on electroweak Su- 
dakov logarithms can be found in Ref. [8]. Specifically, the electroweak corrections 
to the pt distribution of photons and Z bosons at hadron colliders were studied in 
Refs. [HIE]- In Ref. [6], it was found that at transverse momenta of (9(1 TeV) the 
dominant two-loop contributions to these reactions amount to several percent and 
must be included to match the precision of the LHC experiments. 

In this paper we study the electroweak corrections to the hadronic production 
of W bosons at large px- m contrast to the case of Z and 7 production, the 
contributions from virtual and real photons cannot be separated from the purely 
weak corrections and will thus be included in our analysis. 

The partonic reactions qq' — > W ± g (7), q'g — > W ± q (7) and qg — > W ± q l (7) with 
q = u,d, s, c, b are considered. All of them are, however, trivially related by CP- and 
crossing-symmetry relations. Quark- mass effects are neglected throughout, which 
allows to incorporate the effect of quark mixing through a simple redefinition of par- 
ton distribution functions (see Sect. EJ). The calculation of the virtual corrections is 
described in Sect. El In this Section we also present compact analytic expressions 
for the high-energy behaviour of the corrections which include quadratic and linear 
logarithms at one loop as well as quartic and triple logarithms at two loops. The cal- 
culation of the real corrections is performed using the dipole subtraction formalism 
[§1 HOI UI] , as described in Sect. HI After convolution with parton distribution func- 
tions, we obtain radiatively corrected predictions for p T distributions of W bosons 
at the LHC, presented in Sect. El Concerning perturbative QCD, our predictions are 
based on the lowest order. To obtain realistic absolute cross sections, higher-order 
QCD corrections [12] must be included. However, the relative rates for W + , W~ 
and Z production are expected to be more stable against QCD effects. Therefore, 
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in Sect. |5]we also study the impact of the electroweak corrections on these ratios. 



2 Lowest order and kinematics 

The pt distribution of W bosons in the reactions h\h2 — > W ± j is given by 

—, = J2 dx i / dx 2 0(xix 2 - r min )/ /lli i(xi,/i 2 )/ /l2J (x 2 ,/i 2 ) — . 

dp T ~Jj.Ja Jo -■ dp T 



where f min = (p T + m T ) 2 /s, m T = Jp\ + and y/s is the collider energy. The 
indices i,j denote initial-state partons and fh u i(x, /i 2 ), fh 2 j(x,fi 2 ) are the corre- 
sponding parton distribution functions, the partonic cross section for the 
subprocess ij — > W ± k and the sum runs over all k combinations corresponding 
to the subprocesses 

u m d n -> W~g, d n u m -> -> W~« m , d n g -> 1^(1™, 

«m# ~> W 7 "^, fifM m -> W~<2 n , (2) 

for production, and similarly for W + production. The dependence of the par- 
tonic cross sections on the flavour indices m, n amounts to an overall factor \V V w I 2 . 
This factor can be easily absorbed in the parton distribution functions of down-type 
quarks by redefining 

3 3 

fh,dm y I^Wm^il fh,dni fh,d m y I^Wmdnl fh,d n - (3) 

n=l ra=l 

The partonic cross sections can then be computed using a trivial CKM matrix 5 mn . 
The Mandelstam variables for the subprocess ij — > W^k are defined in the standard 
way 

s = (Pi + Pj) 2 , i = (pi - p w ) 2 , u = (pj - p w ) 2 . (4) 

Momentum conservation implies s + t + u = Myy, and the invariants are related to 
the collider energy s and the transverse momentum p^ by p\ = tu/s with s = X1X2S. 
The pt distribution for the unpolarized partonic subprocess ij — ► W^fc reads 



dpx 87riVi,-s|t — u 



(5) 



where J2 — j J2 po i Scoi involves the sum over polarization and color as well as the 
factor 1/4 for averaging over initial-state polarization. The factor l/iV^ accounts 
for the initial-state colour average. 

The unpolarized squared matrix elements for all partonic processes relevant for 
W + and W~ production are related by crossing- and CP-symmetry relations. Thus 
the explicit computation of the unpolarized squared matrix element needs to be 
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performed only once, e.g. for ud — > W g. For this reaction, to lowest order in a 
and as, we have 



ud— >W~g\2 
Born 



32tt z — a s 



a i 2 + u 2 + 2M^s 



til 



(6) 



where s x 



1 — cL denotes the sine of the weak mixing angle. 



3 Virtual corrections 



The one-loop diagrams were reduced to a minimal set of coupling structures, stan- 
dard matrix elements and scalar integrals as in Ref. [6]. The electroweak coupling 
constants were renormalized in the G^-scheme, where a = y/^G^M^s^/n is ex- 
pressed in terms of the Fermi constant G M and = 1 — M^/Mf . Soft and collinear 
singularities resulting from virtual photons were regularized and combined with cor- 
responding singularities from real photons as described in Sect. HI Complete analytic 
results for the one-loop corrections and their asymptotic behaviour will be provided 
in Ref. [13J. The numerical evaluation and detailed cross checks were performed 
with two independent programs. For the scalar loop integrals we used the Fortran 
library [J4] and the FF library |15j . 

In the following, we present compact analytic expressions for the one- and two- 
loop NLL contributions at high energy. As in the case of Z and 7 production [6], 
the NLL terms are obtained from the Born result by multiplication with a global 
factor. For the process ud — > W~g we have 



ud^W~g\2 



Born 



1 + 



2ttJ 



+ 



a 
2^ 



(7) 



At one loop, the NLL part consists of double- and single-logarithmic terms and reads 

I (1) NLL X~YOW fx 2 



3L. 



2 s 2 



U - H 



Here C™ 



loops we obtain 



y q l/(^ 2 w ) + Cf/s 2 w , C f = 3/4, C A = 2, and L 



. ( 2 ) NLL }_f c 

2 V ' 
1 

+ 6 



ew 
6l 



6L?) + 



log(|f|/A4- 

C*A / t 4 , T 4 T 4 



(8) 
At two 



2 s 2 



(14+14 



(9) 



include the 



where 61 = -41/(6(4) and fo 2 = 19/(6s^). The NLL results ©- 
full electroweak corrections in the Mz = Mw approximation. The photonic sin- 
gularities are separated using the fictitious photon mas d3 A = M w . Eq. (ED has 



1 For a discussion of this prescription, as well as for details concerning the implementation of 
the angular-dependent part of the NLL terms, we refer to Ref. [5]. 
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been derived from the general results for leading- and next-to-leading electroweak 
two- loop logarithms in Ref. [3]. 

4 Real radiation 

Soft and collinear singularities from real radiation are treated using the dipole sub- 
traction formalism (9j [TQJ, [TT]. Within this framework infrared singularities of a 
squared amplitude for real radiation are subtracted by an auxiliary function that 
has the same singular behaviour. This ensures that the phase space integral of the 
difference is a finite quantity and the integration can be performed numerically. The 
analytical result for the integral of the auxiliary function over the subspace of a ra- 
diated particle is then added to the result for virtual corrections. The singularities 
of the virtual corrections cancel against those of the integrated subtraction part. 

The algorithms for constructing the auxiliary subtraction function and its in- 
tegrated counterpart have been developed both for the case of photon radiation 
off massless or massive fermions [3] and QCD radiation off massless pU] or mas- 
sive partons [TT]. The latter [TO], [H] approach can be easily adopted to the case 
of photon radiation. It employs dimensional regularization to regularize soft and 
collinear singularities while the former approach [U] introduces small photon and 
fermion masses. We have performed independent calculations of the virtual and 
real corrections within the two regularization schemes and verified that the soft and 
collinear (apart from initial-state and final-state collinear) singularities cancel. We 
also checked that the numerical results obtained within the massive and massless 
regularization scheme are in agreement. In both approaches we use expressions for 
the emission off a massive fermion to describe the emission off a W boson, since 
only soft singularities are present in this case and they depend only on the charge 
of the external particle and not on its spin. 

We restrict our analysis to M^-boson production at high p-r by cutting away events 
with low Pt{W). Additionally, in order to avoid soft-gluon singularities which can 
arise as a side-effect of hard-photon radiation, we introduce a cut on px(jet). 

The remaining initial-state collinear singularities are absorbed into the definition 
of parton distribution functions (PDFs). We choose to perform calculations in the 
MS factorization scheme with the scale M^. The photon- induced processes are not 
included in our calculations since they are expected to be suppressed due to the size 
of the electromagnetic coupling^]. 

2 In order to consistently include 0(a) corrections in a calculation of a hadronic cross section, 
PDFs that are used in the calculation need to take into account QED effects. Such PDF analysis has 
been performed in [TB] and the 0(a) effects are known to be small, both concerning the change in 
the quark distribution functions (below 0(1%) [TT]) and the size of the photon distribution function. 
These effects are below typical uncertainty of hadronic processes. Moreover, the currently available 
PDFs incorporating <D{a) corrections, MRST2004QED [16], include QCD effects at the NLO in 
as- Since our calculations are of the lowest order in QCD, and QED effects on PDFs are estimated 
to be small, we choose to use a LO QCD PDF set without QED corrections incorporated, rather 
than MRST2004QED. 
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For the gluon-induced subprocesses, final-state collinear singularities are also 
present. These singularities from collinear photon-quark configurations are avoided 
by recombining photons and quarks in the collinear region. In practice, we define 
the separation variable 

R(q, 7) = vH-^ + ^-^t) 2 , (10) 

where rji is the pseudo-rapidity and <pi is the azimuthal angle of a particle i. If 
R{Qi 7) < Rsep, then the photon and quark momenta are recombined into an effective 
momentum p T (jet) = Pt(q + 7) which is subjected to the aforementioned cut. 

5 Predictions for the hadronic W production at 
high transverse momentum 

For the numerical evaluation of the corrections we use the following input pa- 
rameters: G M = 1.16637 x l(T 5 GeV~ 2 , M w = 80.39 GeV, M z = 91.19 GeV, 
M t = 171.4 GeV, Mh = 120 GeV. Light-fermion and 6-quark masses are neglected. 
The hadronic cross sections are obtained using LO MRST2001 PDFs [18J. We choose 
fi 2 = p\{W) as the factorization scale and, similarly, as the scale at which the run- 
ning strong coupling constant is evaluated. We also adopt, in agreement with the 
value used in the PDF analysis, the value as(M^) = 0.13 and use the one-loop run- 
ning expression for a>s([i 2 )- For the numerical values of elements in the CKM quark 
mixing matrix we refer to [19J. The statistical accuracy of the Wj cross section at 
LHC is estimated using the integrated luminosity C = SOOftT 1 , and the branching 
ratio BR(W eu e + /i^)=2/9. 

We apply the cuts p T (W) > 100 GeV and p T (jet) > 100 GeV. The value of the 
separation parameter below which the recombination procedure is applied is taken 
to be R sep = 0.1. 

The transverse momentum distributions in LO approximation for pp — > W + j and 
pp — > W~j at y/s = 14TeV are shown in Fig. QJi. In Fig. [Tb we plot the full 0(a) 
electroweak (NLO), one-loop next-to-leading logarithmic (NLL) and next-to- next- 
to-leading order (NNLO) corrections for W + production. The NNLO corrections are 
defined as the NLO plus the two-loop NLL contributions ([9]). As expected, the NLO 
corrections increase significantly with p^. They result in a negative contribution 
ranging from —15% at p^ = 500 GeV to —42% at p^ = 2TeV. The one-loop NLL 
approximation ([8]) is in good agreement (at the 1-2% level) with the full NLO result. 
The difference between NLO and NNLO curves is significant. As can be seen from 
the plot the two-loop terms are positive and amount to +3% at p^ = 1 TeV and 
+9% at px = 2 TeV. The behaviour of the relative corrections for W~ production 
(Fig. [Tb) is qualitatively and quantitatively very similar. 

To underline the relevance of the large electroweak corrections at the LHC, in 
Fig. [2h and Fig. [2b we present the relative NLO and NNLO corrections to the W + 
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Figure 1: Transverse momentum distributions for P^-boson production at the LHC: 
LO predictions (a) and relative electroweak corrections for W + (b) and W~ (c) 
production in NLO (dotted), one- loop NLL (thin solid) and NNLO (thick solid) 
approximation. 
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Figure 2: Unpolarized integrated cross section as a function of p^iW) for W + 
(a) and W~ (b) production: estimated statistical error (shaded area) and relative 
electroweak corrections in NLO (dotted) and NNLO (solid) approximation. 



and W~ cross sections integrated over pr starting from — P™ 1 , as a function of 
p T ut . This is compared with the statistical error, estimated as Acr stat /er = 1/y/N 
with N = £ x BR x er LO where the BR accounts for the ez/ e and [iv^ decay modes 
(for this estimate we ignore experimental efficiencies and cuts). It is clear that 
the size of the NLO corrections is much bigger then the statistical error. Also the 
difference between NNLO and NLO corrections, due to two-loop logarithmic effects, 
is significant. In terms of the estimated statistical error, these two-loop contributions 
amount to 1-3 standard deviations for p^ = 0(1 TeV). 

Finally let us discuss the ratios of the p^ distributions for W + , W~ and Z 
bosons [6]. In contrast to the distributions themselves, these ratios are expected to 
be relatively insensitive to QCD corrections and theoretical uncertainties associated 
with «g and PDFs. Therefore they are good candidates for precision measurements. 
The W + /W~ ratio is presented in Fig. [SK In the considered range the LO value 
increases from 1.5 to 3.5. As already observed, the (relative) electroweak corrections 
to the W + and W~ production processes are almost identical. As a consequence 
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Figure 3: Ratio of the transverse momentum distributions for the processes (a) 
pp^W + j and pp—>W~j and (b) pp—>W + j and pp^Zj at */s = 14 TeV: LO (thin 
solid), NLO (dotted) and NNLO (thick solid) predictions. 

the LO, NLO and NNLO curves in Fig. [3^ overlap. In contrast, the impact of the 
electroweak corrections on the W + /Z ratio (Fig. [3b) is clearly visible. Here the LO 
prediction, ranging from 1.5 to 2, receives corrections that grow with p^ and amount 
to 5-10% for p T > 1 TeV. 

6 Summary 

In this work we evaluated the electroweak corrections to large transverse momentum 
production of W bosons at the LHC. The contributions from real and virtual photons 
cannot be separated in a gauge invariant manner from purely weak corrections and 
were thus included in our analysis. Soft and collinear singularities were regulated 
by introducing a small quark mass and a small photon mass and, alternatively, by 
using dimensional regularization. The real photon radiation was evaluated using the 
dipole formalism. The 0(a) corrections lead to a reduction of the cross section by 
about 15% at transverse momenta of 500 GeV and reach more than —40% at 2 TeV. 
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We also derived a compact (NLL) approximation which includes the quadratic and 
linear logarithms and which gives a good description of the full 0(a) result with 
an accuracy of about 1-2%. Considering the large event rate at the LHC, leading 
to a fairly good statistical precision even at transverse momenta up to 2 TeV, we 
evaluated also the dominant (NLL) two-loop terms. In the high-pp region, these two- 
loop effects increase the cross section by 5-10% and thus become of importance in 
precision studies. We also studied the relative rates for W + , W~ and Z production, 
which are expected to be stable with respect to QCD effects. The electroweak 
corrections cancel almost completely in the W + /W~ ratio. Instead, their impact on 
the W + /Z ratio is significant and amounts to several percent for p T > 1 TeV. 
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